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According to Bell's theorem, the degree of correlation be- 
tween spatially separated measurements on a quantum system 
is limited by certain inequalities if one assumes the condition 
of locality. Quantum mechanics predicts that this limit can 
be exceeded, making it nonlocal. We analyse the effect of an 
environment modelled by a fluctuating magnetic field on the 
quantum correlations in an EPR singlet as seen in the Bell 
inequality. We show that in an EPR setup, the system goes 
from the usual 'violation' of Bell inequality to a 'non- violation' 
for times larger than a characteristic time scale which is re- 
lated to the parameters of the fluctuating field. We also look 
at these inequalities as a function of the spatial separation 
between the EPR pair. 

PACS No(s): 03.65.Bz 



I. INTRODUCTION 

The basic formalism of quantum theory which involves 
concepts of probability amplitudes and the superposi- 
tion principle has been the source of many paradoxes 
and strange implications. A celebrated paradox that il- 
lustrates the counterintuitive conceptual framework of 
quantum mechanics is the 'E-P-R paradox' which is 
based on a gedankenexperiment proposed by Einstein, 
Podolsky and Rosen ( E-P-R) in 1935 [1]. 

E-P-R showed that quantum theory is incomplete as 
the quantum description of the physical system does not 
contain all the 'elements of reality'. The so-called 'ele- 
ments of reality' are those physical attributes of the sys- 
tem that can be measured without disturbing the system. 
E-P-R's argument also required that the principle of lo- 
cality be obeyed in the sense that the disturbance caused 
by the act of measurement does not travel faster than 
light. It seems intuitive to expect all physical theories 
to be local and consistent with realism as stated above. 
But E-P-R's argument showed that quantum mechanics 
is inconsistent with the doctrines of realism and locality. 
A natural requirement at that time was, thus, to look 
for a theory that agreed with all the predictions of quan- 
tum mechanics and yet did not share its conflicts with 
realism and locality. One possibility was to reinterpret 
quantum mechanics in terms of a statistical account of 
an underlying hidden-variable theory which would con- 



form to locality and realism. Many attempts were made 
to find such a complete 'hidden variable theory'. In 1965 
Bell proved that nonlocal correlations among measured 
quantities must obey some general inequalities that are 
common to all local realistic theories [2]. Quantum me- 
chanics, on the other hand, gives rise to correlations that 
violate Bell's inequalities. This provided definite experi- 
mental ways to test for any hidden- variable theories that 
may underlie quantum mechanics. A number of experi- 
ments have been carried out to measure such correlations 
and most of them have yielded results which are in excel- 
lent agreement with quantum theory and in disagreement 
with local realistic theories [3,4] 

In the following we describe Bohm's version [5,6] of 
the E-P-R gedankenexperiment which was considered by 
Bell to prove his theorem. A bound state of two spin-1/2 
particles is prepared in the quantum-mechanical singlet 
state | ijj) given by 
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where \\i n ) and \U n ), i = 1,2, describe the spin states in 
which particle i has spin 'up' and 'down', respectively, 
along the direction n. The bound state is now broken in 
such a way that the particles fly apart in opposite direc- 
tions, but their total spin state remains a singlet. Two 
detectors, far apart from each other , are arranged to 
measure the spin component of particle 1 along a direc- 
tion a and that of particle 2 along a direction b. The 
directions a and b are at an angle 9 with respect to each 
other. According to the quantum formalism, the mea- 
sured expectation value can be written as 



E(a,b) = {ip | cri.ao-2-b | tfj) 



-cos8. 
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One can see that for 9 = 0, E( a , b) = -1. This implies 
the existence of a perfect negative correlation between the 
spins of 1 and 2. Thus, one can predict with certainty the 
result of a measurement on 2 by previously obtaining the 
result of 1. This measurement on 1 does not disturb the 
state of the particle 2, which is very far from 1. There- 
fore, one argues that all components of the spin of 2 are 
'elements of reality', which quantum mechanics cannot 
describe. Further, since the quantum mechanical state 
| ip) does not determine the result of an individual mea- 
surement, there must exist a more complete specification 
of the state through a hidden- variable theory. By assum- 
ing a suitable condition of locality, Bell [2] showed that 
for any local, realistic, deterministic theory, the expecta- 
tion value (2) must satisfy a simple inequality, one form 
of which is 
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E(a, b) - E(a, c) |< 1 + E(b,c) 



(3) 



II. THE MODEL 



where a , b and c are three directions along which the de- 
tectors measure the spin components of 1 and 2. One can 
easily see that for the condition 6( a , b) = 0(b, c) = tt /3 
and 6{ a , c) = 2tt /3, quantum mechanics clearly violates 
this inequality, a result verified by many successful exper- 
iments [4, 6-7]. These results establish nonlocality as an 
inescapable 'fact of life' and cast doubts on the notion of 
'elements of reality'. There are many exciting implica- 
tions of quantum nonlocality like quantum cryptography 
and quantum computers. Particularly fascinating is the 
idea of 'teleportation' of a quantum state using E-P-R 
pairs [8]. 

However, we know from everyday experience that non- 
local correlation of the E-P-R kind are not seen in the 
macroscopic physical world. Classical systems are known 
to conform to realism and locality. How does this transi- 
tion take place? How do nonlocal correlations disappear 
and how does classicality emerge? The transition from 
quantum to classical is marked by the loss of quantum co- 
herence which reduces a pure state density matrix for the 
system to a classically interpretable diagonal mixture. In 
recent years this 'decoherence' [9,10] approach has pro- 
vided a way to understand the emergence of classical- 
ity from underlying quantum dynamics. The approach 
recognizes that a macroscopic system is never isolated 
from its environment but is constantly interacting with 
it. In this approach one calculates the combined den- 
sity matrix of the system and the environment but since 
one is interested only in monitoring the system's degrees 
of freedom, one traces over the environmental degrees of 
freedom. Then one finds that the reduced density matrix 
of the system is driven to a diagonal form in time. 

In this paper we study the effect of an environment on 
an E-P-R singlet where the particles are acted upon by 
forces which push them apart in opposite directions. The 
effect of the environment here is two- fold. Firstly, the two 
particles are coupled to a bath of harmonic oscillators via 
a Caldeira-Leggett type of dissipative coupling through 
their positional degrees of freedom [9] . This aspect of the 
model has been discussed in detail elsewhere [10]. Sec- 
ondly, the spin degrees of freedom for the two particles 
are coupled to an environment modelled by a fluctuating 
external magnetic field. We calculate the spin correla- 
tions for this system, and we show that the system goes 
from a 'violation' of the inequality at t=0 to a 'nonvio- 
lation' over a time scale related to the parameters of the 
fluctuating magnetic field. We also look at this transi- 
tion as a function of the spatial separation between the 
E-P-R pair. In Section II we introduce the model and set 
up the equations that the reduced density matrix of the 
system obeys. In Section III we obtain the solutions for 
these equations and calculate the correlations involved in 
the Bell inequality Finally, in Section IV we summarize 
the results of this paper. 



Our model consists of two particles whose initial state 
can be described by the following wave function 

V = 0(1)^(2) ® -t= (|Tl„>® U 2n >- Uln>® IT2n>) , (4) 

where 0(1 ), <p(2) represent the spatial parts of the wave- 
function which are initially Gaussian wave packets with 
zero initial momenta and centered around x=0 
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The particles are coupled to an environment consisting 
of a collection of harmonic oscillators via their position 
coordinates. Such a model of the environment has been 
studied in great detail by many authors in the context of 
quantum dissipative systems and the quantum measure- 
ment problem [9-10]. In addition to this coupling, the 
spin degrees of freedom for the two particles are coupled 
to an environment modelled by a fluctuating magnetic 
field. The Hamiltonian for this model is: 



p2 p2 

2m 2m 

\ - pj k m lk w 2 lk 
•^2^7 ' —2— [ - ; 

+ E^ + ^l-2 fe 



2m 2 k 



mikLoj k 
C2kX 2 V 



1 



--gh(B 1 .a 1 + S 2 .cr 2 ), 



(6) 



where Pi , X\ and P 2 , X 2 are the momentum and posi- 
tion coordinates of particles 1 and 2, respectively, e is 
the strength of the applied force (the force being positive 
for 1 and negative for 2 ), p\k,x\k and p2k,X2k are the 
momentum and position coordinates of the k th harmonic 
oscillator of the baths which couple to 1 and 2, c\k and c 2 fc 
are the respective coupling strengths and wife and u)2k are 
the frequencies of the oscillators. B\ and B 2 are the ex- 
ternal stochastic magnetic fields acting on particles 1 and 
2, and g is the gyromagnetic ratio. Though in our model, 
we assume that the external fluctuating fields B x (t) and 
B 2 (t) are random in time, we envisage the physical origin 
of this randomness to arise due to spatial variation of the 
field in the region where the singlet is separating. Since 
the analytical treatment for the spatially varying field is 
very difficult, we model the spatial variation seen by the 
wavepackets of the particles to be a temporal variation. 
Bi(t ), B 2 (t) are taken to be Gaussian random processes, 
having the following correlations: 

g 2 {B ia (t 1 )B i!3 {t2)) = I S afj S(h - 1 2 ), (7) 

where i — 1,2 and (3 refer to the components of the vec- 
tor Bi , T a 

oc are the correlation times.. The reduced den- 
sity matrix of the system, p(x\, s\, y\, s^, X2, s 2 , y 2 , s 2 ) = 
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{x\,s\;x 2 ,s 2 | Pr | j/i, Si; J/2, $2) can be obtained by us- 
ing the Feynman- Vernon functional integral method in 
which one first writes the path-integral expression for 
the complete density matrix and then traces over the 
oscillator degrees of freedom. Under conditions of weak 
coupling to the oscillator variables, (i.e., the Markovian 
limit) one finds that pr obeys an equation which can be 
written as 



dpR 
dt 



-i[L sl + L s2 ] P R + ig/2[B 1 .a + B 2 .cr, p R }. (8) 



Pr(xi, si, 2/1, Si; x 2 , s 2 ,y 2 , s 2 ,t) = -p Sl (xi,yi,t) 

Ps 2 (x 2 ,y 2 ,t) 

x[l + CTl .Pi(i)] Slsi 
[l + a 2 .P 2 (t)} S2S , 2 (15) 

In the next Section we write down the solutions for equa- 
tions (12) and (14) and calculate the correlations involved 
in the Bell inequality using Eq. (15). 



Here L Sl and L S2 are a form of quantum-Liouville oper- 
ators which in spatial representation are given by 



III. THE DENSITY-MATRIX OF THE E-P-R 
PAIR 
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where the upper sign is for particle 1 and the lower one 
for particle 2, and 7 and D are related to the environ- 
ment, playing the role analogous to friction and diffusion 
coefficients, respectively. Since the particles are nonin- 
teracting (apart from the correlation due to the initial 
condition) the density matrix equation is separable, i.e., 
Pr = pi ® p 2 with each of p\ and p 2 satisfying the equa- 
tion 

?£i = -iL sjPj + | n r -,-p, -j = 1,2. (10) 

We can further effect the separation between the space 
and the spin degrees of freedom for each particle by rec- 
ognizing that pj is a 2x2 matrix in spin indices and can 
therefore be written as 



The solution for the spatial part of the density matrix, 
i.e., the solution of Eq. (12) have been worked out in 
detail elsewhere [10] . For the initial condition of Eq. 
(5), these solutions are 
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where the upper signs are for j —I and the lower signs 
are for j=2, r = jt, Rj = (xj + yj )/2, rj = Xj — yj, and 



p3 = \[ps J (-v r ii l .in - WjUj.uj.n.rTj 



(11) 



where 1 denotes a 2x2 unit matrix. Substituting Eq. 
(11) in Eq. (10) one easily finds the dynamical equations 
governing p S j and Wj : 



dt 



dt 



iL S jWj — gBjxWj. 



(12) 



(13) 



From Eq.(ll) it is clear that Wj(xjyj;t) = p S j (x 
j,yj,t)Pj and the polarization vector P obeys the equa- 
tion 



dt 



= aPjxBj 



(14) 



Thus the complete solution for the problem can be writ- 
ten as 



M (t) = d 2 + 
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We first note that the off-diagonal part of p S j correspond- 
ing to rj vanishes rapidly in time. If we look at the 
diagonal parts of (16), i.e., the position distribution func- 
tions in the t — > 00 limit, these solutions represent two 
Gaussian wave packets centered around +7^7 and —-^1 
which are moving away from each other with time. 

We now consider Eq. (14). Note that it is exactly the 
equation for a classical moment in an external magnetic 
field. Since the field B(t) is stochastic, P(t) will also 
be a stochastic variable. The quantity of interest corre- 
sponding to a real physical observation would, therefore, 
be the expectation value of P(t), averaged over the en- 
semble of the random process B(t). One has to, thus, 
find the solution to a stochastic Liouvillc equation , i.e., 
an equation for the probability distribution f(P,t) for 
the stochastic variable P{t) [11] . This problem has been 
studied extensively in the literature [11,12] . For the sake 
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of completeness, we record the main results. The equa- 
tion obeyed by the probability distribution function / (P, 
t) averaged over the ensemble of the random process is 



df 
dt 



-Gf, 



(18) 



where 

G = l/2[L 2 x /2r xx + L 2 y /2r yy + L 2 z /2t zz }, (19) 
and the operator L is given by 

d 



-iPx 



dp 



(20) 



From Eq. (18) one can write down the equation for the 
average value of P as 



9(P) 
dt 



= -(G + P), 



(21) 



which reduces to the well-known Bloch equations if wc 
choose t xx = T yy = t\ and t zz — To where the z-axix is 
chosen to be the axis of motion of the E-P-R pair. Eq 
(21) then yields 



d(P z ) 
dt 
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where P± — P x ± iP y . The solutions to the Bloch equa- 
tions (22) are 

(P z (t)) = {P z {0))exp{-t/2 Tl ) 1 (23) 



(P±(t)) = (P±(0))exp(-t/A n - t/AT ), (24) 

Let us now get back to the problem of the E-P-R singlet 
state (1) where the two spins are separately interacting 
with the external fluctuating magnetic fields. We note 
that the spin part of the density matrix corresponding 
to the initial singlet state can be written in the following 
manner: 



P (0) = | V)(V> I 

2\\0 OJ^O 1 

l\ f 
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which can be expressed in terms of the Pauli spin oper- 
ators of particles 1 and 2 as 



p(0) = 1112,1- afa% - u\a\ - crfo-f 



= - ( 1 ® 1 - <7l.<72 



(26) 



This is a convenient form for the initial condition, as we 
note that for calculating the correlation between particles 
1 and 2, the spin part of our solution in Eq. (15) can be 
written as 



Pspin — ^ [ 1 1 



l+<7l.Pl(i)CT2.P 2 (t) 
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Thus the initial conditions can be regarded as the su- 
perposition of three initial condition corresponding to 
Pi = x = -P 2 ;Pi = y = -Pi and P x = z = -P 2 . 
Thus, solving for the time evolution of p involves solving 
for the various components of P{ as in (23), (24)) for 
each of the spins with these initial conditions separately 
and substituting into Eq.(26). One can easily see that 
the time dependent solution is: 

p(t) = Ul ® 1 - Kof + alaDexpi-t/Ar^ 



-t/^-t/ArP-t/^) 
-^a|exp(-t/2r (1) -i/2^ 2) ; 



(28) 



where Tq~\ Tq 2 \ t[ 2 ^ are the characteristic time 
scales for the particles 1 and 2. We can now calculate 
the expectation values required in the Bell inequalities. 
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(i) 



{a1.n1a2.n2) = -(n lx n 2x +nij / n 2y )exp(--(l/r i 

+l/r 1 (1) + 1/rf)) + l/r[ 2) n lz n 2z 
exp(-t/2r ( $ 1) - t/2T {2) ). (29) 

For simplicity, if we assume that all the time scales are 



the same, i.e., Tq 



(i) 



,(2) 
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r s , then one 



can easily see that the expectation value (29) becomes: 

(ai.nia 2 .n 2 ) = Trace(p(t)ai.nia 2 .n 2 ) 

= —cosdexp(—t/T s ), (30) 

transforming the Bell inequality (3) to 

I E(a, b) - E(a, c) | exp(-t/r s ) < 1 + E(b, c)exp{-t / V s ) . 

(31) 

This shows that as soon as the left-hand side becomes 
smaller than unity, the inequalities are satisfied. For the 
choice of angles 9( a , b) — 9(b, c) — it /3 and 9{ a , c) = 2it 
/3, the system goes from a violation of the inequality to 
a nonviolation at time t = —r s ln(2/3). Using the so- 
lution of the spatial part, i.e., Eq (16), this means that 
the quantum correlations are lost when the separation 
between the detectors is of the order of ^T s ln(3/2). As 
mentioned earlier, the origin of the randomness of the 
magnetic field is being attributed to the spatial random- 
ness of the field as seen by the traversing wave-packets 
of each particle. It is clear that for E-P-R correlations to 
persist for longer and longer times, one requires isotropy 
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or absence of field fluctuations in larger and larger por- 
tions of space. Indeed, it is important to keep in mind 
that this instantaneous quantum propagation of informa- 
tion from one particle to another is contingent upon the 
requirement that the intervening space is free from all 
field fluctuations — a requirement that has to be locally 
ensured. Another odd feature worth mentioning is that 
the decoherence in the spatial part and the spin part are 
completely unrelated. Thus, even when the probability 
distributions of position of the two particles have become 
classical and well-separated, the spin wave-function may 
continue to have nonlocal correlations - a situation one 
would not expect classically. 
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IV. SUMMARY 



To summarize, in this paper we have studied the effect 
of the environment on an E-P-R setup. To simulate an 
E-P-R setup, we assume that the two particles of the E- 
P-R singlet are acted upon by forces which move them 
in opposite directions. We find that environmental in- 
fluence causes decoherence in the quantum correlations 
of the initial singlet state leading to a nonviolation of 
Bell's inequalities, thus restoring the classical character 
of the correlations. We have also tried to correlate this 
loss of coherence to the spatial separation between the 
two particles which are additionally coupled to an envi- 
ronment through the Caldeira-Leggett type of dissipative 
coupling. 
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